Closed-loop liquid- vapor equilibrium in a one-component system 
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We report Monte Carlo simulations that show closed-loop liquid-vapor equilibrium in a pure 
substance. As far as we know, this is the first time that such a topology of the phase diagram has 
been found for one-component systems. 

This finding has been achieved on a two-dimensional lattice model for patchy particles that can 
form network fluids. We have considered related models with a slightly different patch distribution 
in the order to understand the features of the distribution of patches on the surface of the particles 
that make possible the presence of the closed-loop liquid-vapor equilibrium, and its relation with 
the phase diagram containing the so-called empty liquids. 

Finally we discuss the likelihood of finding the closed-loop liquid-vapor equilibria on related models 
for three dimensional models of patchy particles in the continuum, and speculate on the possible 
relationship between the mechanism behind the closed-loop liquid vapor equilibrium of our simple 
lattice model and the salt-induced reentrant condensation found in complex systems. 



The recent advances in the chemical synthesis of col- 
loidal particles monodisperse in shape and size have at- 
tracted a lot of theoretical work to understand their col- 
lective behavior. 0, Q. These new particles can be 
thought as 'colloidal molecules' than can organize them- 
selves by self-assembly into useful structures Q. Patchy 
colloids [4j are very interesting because of the possibil- 
ity of tuning the number, interaction, and local ar- 
rangements of their patches. The structure of the self- 
assembled clusters of the patchy colloids can produce 
therefore novel macroscopic behavior. On the other 

hand, primitive models of patchy particles have been used 
to represent the specific anisotropic interactions between 
globular proteins in solution These systems present 
sometimes reentrant condensation induced by variations 
of salt concentration. This behavior is probably asso- 
ciated to effective charge inversion of proteins in the reen- 
trant regime via binding of multivalent counterionsQ. 
The mechanisms involved in the condensation and the 
dependence of the reentrant behavior with the tempera- 
ture are not well understood, but it is likely that specific 
(patchy-like) interactions 0, H| play an important role. 

Given the complexity of both kinds of problems: self- 
assembling of patchy colloids, and phase behavior of pro- 
tein solutions, it seems important the study of simple 
models of patchy particles, in order to understand the 
influence of the numbers of patches and their distribu- 
tion on the surface of the macromolecules (or colloidal 
particles) on the phase behavior and aggregation mech- 
anisms. Hopefully, understanding the phase behavior of 
simple models will be useful to gain insight in the behav- 
ior of those complex systems. 

Recently it has been found that patchy models with 
two types of patches, A, and B in which only A — A, 
and A — B interactions are attractive can exhibit liquid- 
vapor equilibria (LVE) with liquid phases behaving as 
'empty liquids' For these systems the term empty 

liquid [4] refers to the fact that at low temperatures the 



density of the liquid phase, pl on the binodal becomes 
very small at low. Actually pr, approaches zero in the 
limit T — > 0. The empty liquid phenomenology has been 
found for hard sphere particles that carry on the sur- 
face two patches of type A located at opposite poles of 
the sphere, and several patches of type B elsewhere on 
the surface. The potential energy is defined as a sum of 
bonding interactions between pairs of patches belonging 
to different particles. 

Within this interaction scheme, and assuming for sim- 
plicity that the interaction between two particles is zero 
if neither and AB nor an AA bond between them is 
formed, the phase diagram basically depends on the ratio 
between the bonding energies of the two types of possi- 
ble bonds: y — ^abI^aa- From both Monte Carlo Sim- 
ulation and Wertheim's perturbation theory }il-[l2| it is 
known that the regime of empty liquid may occur for 
y < 1/2. Theoretical approximations predict fl2| also a 
lower limit of this ratio, namely yi = 1/3. Simulation 
results are consistent with the existence of such a lower 
bound (yi > 0) for the empty liquid regime. On the 
other hand, for y — 0, the systems at low temperature 
are expected to form linear chains fDij. which, may give 
rise to orientational ordered (nematic) phases at mod- 
erate densities. 3"3| This orientational order can also 
appear as a second phase transition for values of y < 0.5 
where the liquid-vapor is present. The topology of the 



empty liquid phases is well understood. [9Hl2l| At low 
temperature (where the empty liquid regime occurs) and 
not too low density almost all the A-patches are partici- 
pating in a bond, the liquid phases can be then described 
as a network fluid composed by (almost) straight chains 
of spheres linked by sequences of A A bonds. Each of 
the terminal spheres of the chains is connected to an- 
other straight chain, by an AB bond. In the empty liq- 
uid regime the reduction of the temperature implies an 
increase of the typical length of the A A chains, which 
in turn, produces a reduction of the density of the liquid 
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phase. This is due to the fact that A A bonds for y < 1/2. 

At this point, we can ask ourselves to what extent this 
physical picture remains valid if the nature of the col- 
loid surface is changed, and the two A-patches are no 
longer located at opposite poles, but somewhat closer. 
We can define a bending angle, 9ijk for a triplet of par- 
ticles, k} in which j is linked through AA bonds to 
both i, and fc, as: 
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Taking into account that the size of the patches is sup- 
posed to be small in order to preserve the conditions of 
only one bond per patchQ, the mean value of cos Oijk 
will fulfill: ( cos 9ijk) — cos a aa, where a a a is the angle 
between the vectors that connect the center of the sphere 
with the centers of its A patches. Notice that reducing 
the value of olaa from 180° implies an increase the prob- 
ability of forming rings of particles connected by AA- 
bonds, since the system can reduce its energy by forming 
cycles. Moreover, when one of these cycles is formed all 
the ^.-patches of its particles are bonded, and therefore 
these patches cannot bond other chains, which in turn 
can reduce the connectivity of different parts of the sys- 
tem, and eventually inhibit the appearance of LVE. 

Now taking into account that for q^a = 180° and ap- 
propriate values of y one can find the empty liquid regime 
at low temperatures, and that reducing the value of czaa 
might, in principle, destroy the LVE, it seems interest- 
ing to know what happens between the two scenarios. 
In order to investigate this issue, we have found conve- 
nient to analyze a simple lattice model whose results are 
expected to be representative of what can be found for 
models in the continuum, as it has be shown to occur 
in related models lj], 12 [ ■ The computation of phase di- 
agrams is usually much cheaper in terms of CPU time 
in lattice models than in models defined on the contin- 
uum. In spite of their simplicity, lattice models can often 
capture the essential features of the phase diag rams of 
associating molecules and network fluids ll|, ll7H25| . 

We have considered models on the triangular lattice 
(See Fig. [TJ. Sites on the lattice can be either empty or 
occupied by one particle. Each particle carries two sites 
(of type A) that point to two of the nearest neighbor 
(NN) sites. Two particles, i, j; being NN interact with an 
energy U{j — — e if each one points to the other with one 
of its A-patches (AA bond); if only one of the particles 
points to the other with one of its A patches, we have 
Uij = —ye (AB bond, with < y < 1/2); and otherwise 
Uij = 0. 

We have studied three models (See FIG.[T]), which will 
be denominated as o-, m- and p-2A4B models (by bor- 
rowing the prefixes from the organic chemistry nomen- 
clature for the positions of the substituents in aromatic 
cyclic compounds: ortho, meta and para). The three 




FIG. 1: Sketch of the patchy particles in the three lattice 
models: (a) 0-2A4B model: oaa = 60°; (b) m-2A4B model: 
a A A = 120°; (c) p-2A4B model: a A A = 180°. Thin lines 
represent the underlying triangular lattice. Filled circles rep- 
resent particles (occupied sites). Thick segments indicate the 
orientation of A-patches. 



models differ in the arrangements of the patches. In o- 
2A4B we have a A A = 60°, for m-2A4B a A A = 120°, and 
finally p-2A4B model considers aAA = 180°. Notice that 
the latest model is the natural translation of previous 
works for empty fluids into the triangular lattice. In the 
0-2A4B model cycles of AA bonds are easily formed, the 
smallest possible cycle contains just three particles, in 
the m-2A4B model we can find cycles of six AA-bonds, 
whereas in the p-2A4B models no cycles are possible. 
In order to compute the phase diagram we have made 
use of simulation strategies that were described in previ- 
ous papersfll , 1?]; 18 k namely Wang-Landau multicanon- 
ical simulations |llL l2a] to locate critical points and some 
liquid-vapor coexistence points; and Gibbs-Duhem inte- 
gration to trace the binodals. Some special techniques 
have been developed in the simulation procedures to im- 
prove the sampling, including efficient cluster algorithms 
for the P-2A4B model. Technical details will be published 
elsewhere 



27] 



We have have not found evidence of LVE for the o- 
2A4B model. This seems plausible, since at low temper- 
ature and moderate density the particles can easily ar- 
range themselves in triplets that minimize the energy and 
block the development of an extended network of patch- 
patch bonds. Both m-2A4B and p-2A4B models show 
LVE for some range of values in the region y < 1/2. The 
p-2A4B model shows, as expected, a phase diagram (See 
FIG [2]) that fully resembles those of the analogue models 
in the continuum [t| [Io| , and on the square lattice [ll| ; 
i.e. it appears the regime of empty liquid. 

The key result of this paper is, however, the phase 
diagram of the m-2 A4B model (shown in Fig [3]) . For 
y = 0.50 we get a similar behavior in m-2A4B and p- 



3 




FIG. 2: Temperature-density phase diagrams for the p-2A4B 
model, and different values of y (See the legends). Filled sym- 
bols indicate critical points. T* is the reduced temperature: 
T* = fesT/e, and the density p corresponds to the fraction of 
occupied sites. 
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FIG. 3: Temperature-density phase diagrams for m-2A4B 
model. See FIG. [2] caption for notation. 



2A4B models, i.e, we find typical LVE binodals. How- 
ever, in the m-2A4B model we find that the LVE for 
some range of values of y ( yi < y < 1/2), with y\ ~ 0.43, 
exhibits two critical points, that will be named, according 
to their critical temperature as lower critical point (LCP) 
and upper critical point (UCP). The LVE only occurs 
between these two temperatures, which of course depend 
on the value of y. The critical points were determined 
through a finite-size scaling analysis on the multicanon- 
ical simulation similar to that explained in Ref. 



sidering typically system sizes in the range 12 < L < 48 
for the m-2A-4B model (with L 2 being the total number 
of sites) . The scaling behavior was found to be fully con- 
sistent with the two-dimensional Ising criticality for both 
LCP and UCP. 

The existence of reentrant phase equilibria is well 
known to appear in some binary mixtures. However, for 
one-component systems it is very unusual to find fluid- 



fluid reentrant critical behavior. Roberts et al. 20} found 
such a behavior for a lattice model of network-forming 
fluids: However the closed-loop coexistence found in Ref. 
[20l ] involves liquid phases. In our case, the phase equilib- 
ria can be appropriately described as LVE, by considering 
that modifying the value of the interaction parameter y 
the reentrant phase diagrams can be connected continu- 
ously to typical LVE diagrams that appear for y > 1/2. 
Along this path the LVE envelope varies smoothly for 
T > 0. 

In the previous paragraph we have indicated that the 
closed-loop LVE occurs for any value of y in the range 
(yi, 1/2), however our simulation results only support, in 
principle, the reentrance for yi < y < 0.46. As y ap- 
proaches 1 /2 it becomes a very hard task to compute the 
LCP using Monte Carlo simulation, since it occurs at 
very low temperatures. However, we can make use of the 
scaling criteria introduced in Ref] 111: At very low temper- 
ature, the probability of finding non-bonded A patches 
becomes negligible small, therefore we can consider as a 
good approximation to analyze the phase equilibria con- 
sidering only those configurations in which each patch A 
forms either an AA or an AB bond. Under this assump- 
tion the potential energy of the model can be written 
as: [ll1| (U/e) = -N + Nab (5 - y), where N is the num- 
ber of occupied sites and Nab the number of AB bonds, 
Therefore the phase behavior at low temperature, T, will 
basically depend on the ratio (1/2 — y)t/ksT (with ks 
being the Boltzmann constant). Then the temperature 
at the LCP is expected to scale for (1/2 — y) — > + as: 
T l c oc (1/2 — y). The low temperature phase diagram of 
the m-2A4B model has been found to fulfill very well this 
scaling at low temperature (See FIG. HJ). Of course the 
reliability of this prediction depends on the absence of 
additional possible phase transitions that could occur in 
the model. Notice that the p-2A4B model is expected to 
exhibit and order-disorder transition, llj, ll4H16| Such a 
transition will merge with the LV transition as T — > 0. In 
the case of the m-2A4B model we have not found any sig- 
nature of a second phase transition at the temperature 
range where the close- loop LVE occurs. The origin 
of the different topologies of the LVE phase diagrams 
of p- and m-2A4B models can be explained by looking 
at representative configurations of both models at low 
temperature (FIGS I5H7|). In the p-2A4B model the end 
particles of each chain are bonded via ^4i?-bonds to two 
other chains on the system, then on average each chain 

at 



111 con- is bonded to four neighbor chains (Nab — 2iV< 



chains 



low T). At very low temperature in the m-2A4B model 
we can see (See FIG [5]) the formation of rings with differ- 
ent number of particles, and the existence of AB intra- 
chain bonds. As a consequence the number of available 
A patches to form interchain bonds is reduced, and the 
bonds do not produce a macroscopic bond network. No- 
tice that as one increases the temperature (See FIG [7]) 
there is a reduction of the number of rings due to en- 
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FIG. 4: Liquid vapor equilibria phase diagrams for m-2A4B 
model with an alternative temperature scaling. 
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FIG. 6: Configuration for the m-2A4B model with p — 1/3, 
T* = 0.03 (below the lower critical temperature). 




FIG. 5: Configuration for the liquid phase at the binodal for 
the p-2A4B model with y = 0.40, T* = 0.07, p ~ 0.375 
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FIG. 7: Configuration for the m-2A4B model with p = 1/3, 
T* = 0.10. At these conditions there is LV phase separation 



or 



tropic effects, [28| this enhances the connectivity and the 
phase separation can take pface, by a mechanism similar 
to that present in the p-2A4B model. 

It has to be emphasized that the origin of the reen- 
trance of the LVE is not just the angle between patches 
in the colloidal particles, but the fact that at low temper- 
ature the particles have the trend to self-assemble form- 
ing small aggregates (usually rings) in which most of their 
sticky spots are saturated by intra-aggregate interactions. 
This process might produce a reduction of the effective 
valence of the aggregates with respect to that of the inde- 
pendent particles, which in turns reduces the possibility 
of bond network formation and subsequently hinders the 
ability of the system to condense. 

The mechanism underlying the reentrance of the LVE 
is rather different to that appearing in binary liquid mix- 
tures, in which the reentrant solubility at low tempera- 
ture is attributed to the effects of large anisotropic in- 
teractions (hydrogen bonding) between the two compo- 



nents, whereas the demixing at intermediate tempera- 
tures is due to Wan der Waals interactions, which fa- 
vor the interaction between pairs of particles of the same 
species. IS 3(| 

Notice that the results presented in this communica- 
tion resemble, to some extent, those reported by Sciortino 
for a system of Janus particles, and by Rein- 



et al. 31 



hardt et al. 33] for systems with competition between 
phase separation and self-assembly. In both cases, at 
low temperature, the vapor density at LVE increases on 
cooling. However in these systems no LCP appears, due 
to the presence of crystalline phases that make both liq- 
uid and vapor metastable prior the possible presence of 
the LCP. 

The question that arises after having found the closed- 
loop LVE in the lattice model for systems with two A- 
patches with olaa = 120°, is to what extent these results 
can be extrapolated to models in the continuum, either 
in two or three dimensions. In principle, one can expect 
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that for 2D-models the formation of rings, which in turn 
may produce the closed-loop LVE, will be enhanced for 
models with angles cxaa that fit those of regular polygons 
q.aa — I ir^ 7r J w hh n = 4, 5, 6, • • • . For 3D models one 
can also expect the formation of rings, in this case the 
presence of additional degrees of freedom in the chains 
will give the systems more flexibility for the formation 
of rings, however, on the other side, the larger space di- 
mensionality of the system will favor the presence of open 
chains due to entropy effects. In any case the pursue of 
the conditions whether closed -loop LVE will be obviously 
not so straightforward as in the current lattice models. 
We plan to work on these issues in the future. 

As a summary, we have shown how patchy models on 
the triangular lattice can exhibit closed-loop LVE on one- 
component systems. The basic physics underlying this 
behavior is the reduction of valence at low temperature 
by the formation of compact structures at low and mod- 
erate densities that reduce the amount of connectivity 
between different parts of the systems. The phase di- 
agram depends dramatically on the angle between the 
A patches. As far as we know, the work presented here 
constitutes the first model exhibiting closed-loop LVE on 
a one-component system. Finally, self assembling mech- 
anisms similar to those appearing in the simple models 
presented in this paper, might play some role in the reen- 
trant condensation phenomena of globular proteins [1] 
and polyelectrolytes [34j. 
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